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This paper presents dynamical models of a large flexible launch vehicle. A complete set of 
coupled dynamical models of propulsion, aerodynamics, guidance and control, structural 
dynamics, fuel sloshing, and thrust vector control dynamics are described. Such dynamical 
models are used to validate NASA’s SAVANT Simulink-based program which is being used 
for the preliminary flight control systems analysis and design of NASA’s Ares-1 Crew 
Launch Vehicle. SAVANT simulation results for validating the performance and stability of 
an ascent phase autopilot system of Ares-1 are also presented. 


Nomenclature 

R e = Earth’s equatorial radius = 20925646.325459 ft 

R p = Earth’s polar radius = 20855486.595144 ft 

J 2 = Earth’s second order zonal coefficient = 1.08263 lx 1(T 3 

J, = Earth’s third order zonal coefficient= —2.55x10 6 

J 4 = Earth’s fourth order zonal coefficient = —1.61x10 6 

U = Earth’s gravitational potential 

H = Earth’s gravitational parameter = 1 .407644 176x1 0 16 ff / .S' 2 

(j) = Earth’s geocentric latitude 

( Sx j S y j 8- ) = ( X >Y> Z ) components of the gravitational acceleration 

r = vehicle’s position vector 

r = magnitude of vehicle’s position vector 

( X,y,z ) = (x,y,z) components of vehicle’s position vector in an inertial reference frame 
( a x , a y , a z ) = (x,y,z) components of vehicle’s absolute acceleration in an inertial reference frame 

V = vehicle’s absolute velocity vector 
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(u,V, w) = (x,y,z) components of vehicle’s absolute velocity vector in an inertial reference frame 


V, = relative velocity vector of the vehicle (measured in a body-fixed reference frame) 

Q = angular velocity vector of the Earth 

( y = z component of the Earth’s angular velocity vector = 7.2921 152087 x 10 5 rad/s 

CO = angular velocity vector of the vehicle 

(p.qjj = (x,y,z) components of vehicle’s angular velocity vector in a body-fixed reference frame 


V w = velocity of the wind 

V = air-stream velocity vector 

V = magnitude of air-stream velocity 

(V,„ xt’Kn vb’Kn zb) = ( X >Y> Z ) components of air-stream velocity in an inertial reference frame 
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= undamped natural frequency of the actuator dynamics= 32.6726 rad/sec 
= pitch gimbal angle 
= pitch gimbal angle command 
= pitch gimbal angle acceleration 
= yaw gimbal angle 
= yaw gimbal angle command 

= yaw gimbal angle acceleration of the actuator dynamics 

= total thrust inside the atmosphere 
= total vacuum thrust 

= nozzle exit area = 122.137/?" 


p 0 = local atmospheric pressure 
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( F aero.xb > F aero.yb » F aero. zb ) = ( x >y> z ) components of aerodynamic force in Body Frame 
( F rkt. xb ’ F rkt.yb > F rkt.zb ) = ( X ’Y’ Z ) components of rocket engine force in Body Frame 
( F r CS . xb , F rcs.yb > ) = ( X >Y> Z ) components of reaction control force in Body Frame 

( F sbsh.xb » » F ziosh.zb ) = ( X >Y> Z ) components of slosh force in Body Frame 
> F total. z b ) = ( X >Y> Z ) components of total force in Body Frame 
( F ,otai.xi » F total. yi » F to,ai.zi ) = ( X >Y’ Z ) components of total force in Inertial Frame 
(^.*6 , F aero.yb > ) = ( X >Y> Z ) components of aerodynamic torque in Body Frame 

(T rfa ^ , T rkt yb , T rkt zb ) = (x,y,z) components of rocket engine torque in Body Frame 

(7) cs x6 , T rcs yb , T rcs zb ) = (x,y,z) components of reaction control torque in Body Frame 

( T sbsh.xb » T siosh.yb » T slosh. zb ) = ( x >y> z ) components of slosh torque in Body Frame 
T TWDp = pitch torque on the vehicle due to the TWD effect 

T TWDy = Y aw torque on the vehicle due to the TWD effect 

/ = nozzle inertia in plane of movement = 19102.0833 Ibf ■ ft ■ S 2 

M e = nozzle mass = 694.4 lb 

/ = distance from vehicle’s center of gravity to nozzle pivot point 

l = distance from nozzle pivot point to nozzle center of gravity = 1 .2775 ft 

r = vector from vehicle’s center of mass to slosh fuel center of mass in Body Frame 

/ = center of mass position vector of vehicle in Body Frame 

/ tank = tank location vector in Body Frame 

/ = slosh moment arm 

M s = slosh mass 

£ = damping ratio of the slosh fuel dynamics 

CO s = undamped natural frequency of the slosh fuel dynamics 

rj = flex mode state 

i^l lex = damping ratio of flex modes 

OJj lex = undamped natural frequency of flex modes 

171 = vehicle mass 

(c t , C v , C, ) = (x,y,z) components of center of mass 
T s = sampling period = 0.02s 

K p = proportional gain 
K j = integral gain 

K d = derivative gain 

I. Introduction 

N OTE to Session Organizer/Reviewers: This draft manuscript summarizes the preliminary results obtained 
during an early phase of a new project for the dynamical modeling and flight control design of large flexible 
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launch vehicles as applied to Ares-I Crew Launch Vehicle. During the next several months, a more detailed, 
rigorous study will be conducted in the areas of coupled dynamical modeling of propulsion, aerodynamics, guidance 
and control, and vehicle structure. A companion paper on flight control systems analysis and design for large 
flexible launch vehicles is also being submitted to the Space Exploration and Transportation GNC session. 

II. Definition of Coordinate Frame 

A. Geocentric Equatorial Inertial Frame or Inertial Frame 

Geocentric equatorial inertial frame or simply inertial frame (Fig. 1). Origin is at Earth Center. Axis z. is normal 

to equatorial plane, pointing to North Pole; Axes X and y f are in equatorial plane, axis X i is along direction of vernal 
equinox, which is the direction of intersection of Earth equatorial plane and Sun ecliptic plane. 

B. Geocentric Equatorial Rotating Frame or Central Earth Frame or Earth Frame 

Geocentric equatorial rotating frame is fixed to the Earth, also called central Earth Frame (Fig.2). Origin is at 
Earth center. Axis z e is normal to equatorial plane, pointing to North Pole, hence coincides with z ; . . 

Axes X e and _V e are in equatorial plane, with axis X e in Greenwich meridian. This frame has angular velocity of Earth. 


North rotational pole 



Greenwih 

meridian 


North rotational pole 



Figure 1. Inertial Frame 


Figure 2. Earth Frame 
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C. Body Frame or Body-fixed Frame 

Body-fixed frame, simply body frame (Fig.3), is 
rigidly fixed to the vehicle body. Origin is at empty 
vehicle center of mass; Axis X h is along structural 

longitudinal axis, pointing forward; normal axis z, is in 
plane of symmetry, perpendicular to X h and pointing 
downward; axis y h is perpendicular to plane of 
symmetry and pointing rightward. 

III. Aerodynamics Forces and Moments 



From the state variables, we can determine the 
position and velocity of the vehicle in the inertial frame. 

And then, the geographical latitude, the height over the 
Earth surface, angle of attack, angle of sideslip and 
Mach number can be calculated. By looking up the data 
table, aerodynamic coefficients and base force can be 

found. The equations of aerodynamics forces and moments can be written as follows: 
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Figure 3. Body Frame 
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D = C a QS - F base 

< c = c y/3 pqs 

N = (C N0 + C Na a)QS 

Where, C A is axial force coefficient; C Yf! is side force curve slop; C vo is normal force coefficient at zero angle 
of attack; C Na is normal force curve slop. 


F h =-D 

aero.xb 

F -C 

aero.yb 

F h =-N 

aero. zb 
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T a ero,b = c y {C m + C Na a)QS + c z C yp PQS + C Mrp QSb ref 
' Lero.yb = ~ C : ( F base ~ C aQ S ) + {~ARP ~ C X )(C NQ + C Na (X)QS + (C MpQ + C Upa O)QSC ref 
Kem.zb = i-ARP - c x )C yp J 3QS + c y (F base - C a QS) + C myP j3QSb ref 


Where, C Mrp is rolling moment coefficient; C Mp0 is pitching moment coefficient at zero angle of attack; C Mpa is 
pitching moment curve slop; C p is yawing moment curve slop. 

In the SAVANT, C N0 — C Mp0 = C Mrp = 0 . 

IV. TWD Model 

Tail Wag Dog model 

F TWDp = i^e + M e l C gl e )d p 

T TW Dy={I e +MJc g le)S y 

V. Flex Model 

Flex state dynamics 

n + K flex (o aex r\ + 0)^T1 = F r T kt <p rkt 

Where, t] is the flex mode state column vector. F rkt is the transpose of the rocket engine force column vector, 
and (p kt is a 3 by 6 flex mode parameter matrix. 


Sensor error 
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lh2 Hcx _ 
lh2 Hcx _ 9, ,,,2*1 
*Plh2 flex _ <Pslm2n 


Gimbal compliance 


'P iktflcv = Vrk, 1 ! 


Where, \| / navl , \|/ nav2 , \|/ nav3 , <p s|ml , <p s|m2 , \|/ rkt are 3 by 6 parameter matrix respectively. 

VI. Slosh Model 

The sloshing will be modeled as a spring-mass-damper system in y-z plane; we do not consider the x component. 
For this program, it does not include the flex model effect to the slosh model. 

K = - (o)r s - {r rel +ojx (r - I cg - l tmk ) + 2 coxr + d)x[ojx (r - l cg - l tmk )]} 

Matrix Form in Body Frame: 
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VII. Rocket Model 

Rocket Propulsion: 

T = T 0 ~Po A e 

The x, y and z components of the thmst in the Body Frame: 
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The x, y and z components of the torque due to the thrust in the Body Frame: 
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VIII. Gravity model 

r = y/x 2 + y 2 +z 2 
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IX. Actuator Model 

The actuator model is a second order system. 

S + 2 Cm 8 + co~S = (o~8 

p “ n p n p n pc 
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X. Force and Momentum 
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To transfer the (x,y,z) components of total force from Body Frame to Inertial Frame ,by the 
quaternion (t/j q 2 q 3 q 4 ) . 
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In Inertial Frame: 
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Kinematical Equation of Rotation: 
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Angular acceleration: 
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XI. Quaternion 


To renormalize quaternion: 
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Where, (q l q 2 q 3 </ 4 ) is the quaternion column vector and (q in q 2n q 3n q 4n ) is the quaternion 

column vector after renormalization. 

Quaternion derivatives: 
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XII. Guidance Commands 

XIII. Control System 

Discrete PID controller for roll, pitch and yaw channels 

u(t) = K p e{t) + KY i (t) + K d e{t) 
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Discrete integration 


Yj(t) = Y t {t -\) + ^-[e{t) + e{t -\)\ 


1. Center of Pressure 


XIV. System Analysis 

M LE=~c x L + M cg =-x cp L 


Where, M LE is the pitching moment relative to the pin point of the vehicle; M ca is the pitching moment with 
respect to the center of mass, X cp is the location of center of pressure. 

According to the definition (force-and-moment system) of center of pressure, the plot of center of pressure can 
be seen in the following figure. 



1 1 

1 1 


1 ~~ 
CX 


1 1 
1 1 


cpx 

— 

1 1 

t r 

i i 

i i 

1 

1 

1 

— 

" i X i 

i a. 

4 1- 

i i 

i i 

i i 

i i 

i i 

1 

1 

4 

1 

1 

1 

1 

1 


i i 

i i 

i i 

i i 

i i 

t i- _ — 

1 

1 

1 

— 

i i 

i i 

1 

1 

1 

1 


Time, sec 

Figure 4. Location of the center of pressure center of mass 



Figure 5. Vehicle position in Inertia Frame 
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Mach Number 
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Figure 6. Moments of inertia 



Figure 7. Relative velocity 
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Figure 8. Mach number 




Figure 9. Dynamic pressure 
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Angle of Sideslip, deg Angle of Attack, deg Pitch Gj m b a | Angle Command, rad 




Figure 10. Pitch gimbal angle command 


Figure 11. Yaw gimbal angle command 
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Figure 12. 


Angle of attack and sideslip angle 


Figure 13. Total thrust 
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Weight Flow Rate, Ibf/sec 
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Figure. 14 Solid rocket booster weight flow rate 



Figure 15. Solid rocket boost specific impulse 
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Figure 16. Euler angle 
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Figure 17. Axial force coefficient 
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Figure 18. Normal force curve slop Figure 19. Side force curve slop 
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Figure 20. Pitching moment curve slop 


Figure 21. Yawing moment curve slop 
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Total Mass, Ibf 
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Figure 22. Total Mass 


2. Stability Analysis 


XV. Sample Case 


<P rkt 


0.000000272367963 0.000000174392026 -0.000000347086527 
-0.000364943105155 0.006281028219530 0.000491932740239 
0.006281175443849 0.000364891432306 -0.006260333099131 


-0.000000266173427 0.000000329288262 -0.000000369058169" 
-0.006259406451949 -0.000542750533582 -0.007673360355205 
-0.000491798506301 0.007676195145027 -0.000542218216634 y 
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Table 1. Initial States’ Values 


Parameter 

Initial Value (unit) 

U 

1340.649569297657 (ft/sec) 

V 

6.542058392976782 (ft/sec) 

W 

0.244274988189246 (ft/sec) 

P 

0.000034916172322 (rad/sec) 

0 

0.000064018398388 (rad/sec) 

R 

1. 016439536705 16e-20 (rad/sec) 

X 

87898.84842619013 (ft) 

Y 

18384832.09303243 (ft) 

Z 

9960462.310921878 (ft) 

qi 

0.359443352340896 

q2 

0.608859075914155 

q3 

0.362476166996897 

q4 

0.60720847366793 


XVI. Conclusion 
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